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Abstract 

A method for obtaining the low chemical potential dependence of the dressed quark propagator 
from an effective quark-quark interaction model is developed. From this the chemical potential 
dependence of the chiral quark condensate and bag constant are evaluated. It is found that the 
chiral order parameter decreases with the increasing chemical potential jx, which is qualitatively 
different from the results given by some authors in the literature. The reason for this discrepancy 
is discussed. 
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High-energy heavy-ion collision experiments, which produce systems with large baryon 
density, are an important step in the search for quark-gluon plasmafl]. Present efforts 
are concentrated on exploring the domain of nonzero baryon density. The field theoretical 
study of finite baryon density in QCD begins with the inclusion of a chemical potential, 
/i, which modifies the fermion piece of the Euclidean action: 7 • d + m — > 7 • d — 74/i + 
m. The introduction of a chemical potential leads to a domain where the nonperturbative 
information can only be obtained from models(The numerical simulations of lattice-QCD 
action are not tractable in the domain of /j, 7^ 0[2]). Since the global color symmetry 
model(GCM) [3,4,5] provides a nonperturbative framework that admits the simultaneous 
study of dynamical chiral symmetry breaking and confinement, it is expected to be well 
suited to explore the transition from hadronic matter to QGP[6]. 

It is the aim of this letter to study the chemical potential dependence of the dressed- 
quark propagator in the framework of the GCM, which provides a means of determining the 
behavior of the chiral and deconfinement order parameters. Up to this end let us consider 
the Euclidean action of the GCM at finite chemical potential //(In the case of the chiral 
limit): 



the dressed model gluon propagator in GCM. We want to stress here that g 2 s D ab v {x — y) 
does not evolve with //. For convenience, we will employ a model ansatz D"^(x — y) — 
5 fll/ 5 ab D(x—y) for the gluon propagator, which is often referred to as the so-called "Feynman- 
like" gauge[4](It should be noted that the above ansatz should be regarded merely as a model 
form for the gluon two-point function). 

Introducing an auxiliary bilocal field B e (x,y) and applying the standard bosonization 
procedure the partition function of GCM [3,4] 



SGCM[q,q;/4 



= J d 4 x{q(x)fr-d x -wMx)} + J d 4 xd 4 y ^%(x)D%{x -y)j h M , (1) 



where j^(x) = q(x)^^q(x) denotes the color octet vector current and g 2 s D 



y) is 




(2) 



can be rewritten in terms of the bilocal fields B e (x,y) 




(3) 
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with the effective bosonic action 

S eff [B«; A = -Trlng-i[B»; „\ + j ^y ^^f^ (4) 

and the quark operator 

g- 1 [B e ; ^ = [7.9,- n lA ]5{x -y) + A e B e (x, y). (5) 

The matrices A e = D a ®C b ®F c is determined by Fierz transformation in Dirac, color and 
flavor space of the current current interaction in Eq.(l), and are given by 

One might expect that the complete set of the 16 Dirac matrices {Id, 75, 7^, 757^, a^} 
must be employed in the description. However, by limiting the gluon two-point function 
g 2 D ab v [x — y) to diagonal components in Lorentz indices, the tensor a^ v is excluded. 

In the mean-field approximation, the fields B 9 (x,y) are substituted simply by their vac- 
uum value B^(x,y), which is defined as SS ^ S \ B = and is given by 

B e M{x,y) = g 2 D(x - y)tr[A 9 g \pL](x,y)], (7) 

where the notation tr includes trace over the Dirac, color and flavor indices and g$ l (x, y) de- 
notes the inverse propagator with the self-energy S(x, y) = A 9 Bq(x, y) at the finite chemical 
potential /i. Employing the stationary condition Eq.(7), and reversing the Fierz transfor- 
mation, we have 
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£M (2/1,2/2) = 3 g 2 s D(yi,y2hvGo[tA (yi, 2/2)7*- (8) 

It should be noted that both Bq(x, y) and g^^x, y) are dependent on the chemical po- 
tential /i. If the chemical potential /i is switched off, t/oM g° es m to the dressed vacuum 
quark propagator G = (? [/i = 0], which has the decomposition 

G-\p) = i 1 -p + E(p) = i 7 • P A(p 2 ) + B(p 2 ) (9) 

with 

E(p) = / d'xe^\A d Bl(x)} = tj -0^9 2 s D(p - q)-y v G{q) lv (10) 

where the self energy functions A{p 2 ) and B{p 2 ) are determined by the rainbow Dyson- 
Schwinger equation(DSE) 
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Here we want to stress that the B(p 2 ) in Eq.(ll) has two qualitatively distinct solutions. 
The "Nambu-Goldstone" solution, for which 

B(p 2 ) ^ 0, (12) 

describes a phase in which: 1) chiral symmetry is dynamically broken. Because one has a 
nonzero quark mass function; and 2) the dressed quarks are confined, because the propagator 
described by these functions does not have a Lehmann representation. The alternative 
"Wigner" solution, for which 

B(p 2 ) = 0, (13) 

describes a phase in which chiral symmetry is not broken and the dressed- quarks are not 
confined. In "Wigner" phase, the Dyson-Schwinger equation(ll) reduces to: 

where A'(p 2 ) denotes the dressed quark vector self energy function in "Wigner" phase. 
Therefore, the dressed quark propagator in "Wigner" phase can be written as G^ w \q) = 

A'(q 2 )q 2 ■ 

In order to get the numerical solution of A{p 2 ), B{p 2 ) and A'{p 2 ), one often use model 
forms for gluon two-point function as input in Eq.(ll). Here we investigate two different 
two parameter models for gluon propagator; 

^<V) = 9 2 Dfl{q 2 ) + g>D uv (q 2 ) = ^ 2 e-i + v , (14) 

A q 2 ln(j^ + e) 

\ QCD ) 

and 

g 2 D^\q 2 ) = g 2 D%(q 2 ) + g 2 D uv ( q 2 ) = ^d^— + ^ ^ (15) 

\ QCD J 

The term Dj^(q 2 ), which dominates for small q 2 , simulates the infrared enhancement and 
confinement. The dressed-gluon propagator is strongly enhanced, which leads via the QCD 
gap equation(ll) to an infrared enhancement of the light quark mass function. These modifi- 
cation are intimately related to the confinement and dynamical chiral symmetry breaking[7]. 
The other term D uv (q 2 ), which dominates for large q 2 , is an asymptotic ultraviolet (UV) tail 
match the known one-loop renormalization group result with d = [12/(33 — 2Nf)] = 12/27, 



A-qcd = 200 MeV. The model parameters x an d A are adjusted to reproduce the weak 
decay constant in the chiral limit f w = 87 MeV. The forms of g 2 D(q 2 ) have been used in 
Ref. [8] and it has been shown that with these values a satisfactory description of all low 
energy chiral observables can be achieved(more detail can be seen in Refs.[8] and [9]). 

Let us now study the chemical potential dependence of the dressed quark propagator. To 
this, one can numerically study Eq.(8). From Lorentz structure, the most general form for 
the G~ l {p, A*) (in the chiral limit) which fulfills Eq.(8), reads 

G~ 1 {p, /i) = ■ p a(p 2 , u 2 , [u ■ p} 2 ) + b (p 2 , it 2 , [u ■ p} 2 ) + 27 • up ■ u c(p 2 , u 2 , [u ■ p] 2 ) (16) 
-M7 -pp-u d(p 2 , u 2 , [u ■ p} 2 ) — 7 • u e(p 2 , u 2 , [u ■ p} 2 ) +ip ■ u f(p 2 , u 2 , [u ■ p] 2 ), 

where = (0, /x). However, the model gluon propagator(14,15), as pointed out above, has 
no explicit /i-dependence, which can arise through quark vacuum polarisation insertions. 
As such it may be inadequate at large value of /1, particularly near any critical chemical 
potential. Therefore, it is comparatively safe to study the low chemical potential dependence 
of the dressed quark propagator by means of Eq.(8) with the model gluon propagator(14,15). 
Here we restrict ourselves to study the first-order dependence of ^o" 1 ^] upon //. Namely, we 
can expands ^o" 1 ^] in powers of fi as follows 



+ 



which leads to the formal expansion 



fx + 0(fx 2 ) = G- 1 + fxT^ + Oifj 2 ), (17) 



QM = G-GnF 4 G + ---, (18) 

with r 4 



r 4 (yi,y2j = — 



(19) 

fi=0 



5/i 

In coordinate space the dressed vertex T 4 (x, y) is given as the derivative of the inverse quark 
propagator ^o^M with respect to the chemical potential /i. 

Taking the derivative in Eq.(5) and putting it into Eq.(19), we have 

<^Lu](yi,y 2 ) 



r 4 (yi,z/2) = -74^(yi - w) + 



(20) 



5fi 

Substituting Eq.(8) and (18) into Eq.(20), we have the inhomogeneous ladder Bethe 
Salpeter equation(BSE) for T 4 vertex, which reads 

r 4 (2/1,2/2) =-74^(2/1-2/2) 



4 r 

--g s D(y 1 -y 2 ) J du 1 du 2 'y v G(y 1 , UijF^Ui, u 2 )G(u 2 , y 2 )~f v (21) 
Fourier transform of Eq.(21) leads then to the momentum space form of T 4 



r 4 (P, 0) = -74 - \ J 0yi9 2 s D(P - K) lu G(K)T,(K, 0)G(K) lv 



Eqs.(9) and (22) yields 



T 4 (P,0) 



.dG^jP) 
1 <9P 4 : 



(22) 



(23) 



i.e., the vector "Ward identity" is satisfied [10, 11]. This means that it is the nonperturba- 
tive dressed effect which modifies the bare vertex 74 to the nonperturbative dressed vertex 
T 4 (P, 0) at the level of linear response approximation. 

By means of "Ward identity" (23) and Eq.(17), we have the chemical potential dependence 
of the dressed quark propagator in "Nambu-Goldstone" and "Wigner" phase separately ( here 
we only consider the first order dependence of upon /i at the mean field level); 



6f G| M 



«7 • p 
+274 

27 • p 
+il4 



A{P 2 ) + 2ipP 4 



dA{P 2 
dP 2 



P A A{P 2 ) + ^ P f-^p ■ 



+ P(P 2 ) + 2ipP 4 
] -+ipA(P 2 ) 



dB{P 2 ) 



dP 2 
+ 0(p 2 ), 



(24) 



A'(P 2 ) + 2ipP A 



dA'{P 2 
dP 2 



P A A'{P 2 ) + 2ippf-^ ) ' .-"'"^ 



(25) 



Note that the functions A(P 2 ), B(P 2 ) and A'(P 2 ) here, as pointed out above, are known by 
numerically solving the corresponding DSE. 

As it is shown in Eqs.(24) and (25), for p ^ the dressed-quark self energies in general 
acquire an imaginary part driven by the chemical potential /z[12-15]. Just as pointed out 
in Ref.[6], this effect is not observed in the study of the Nambu-Jona-Lasinio model and 
the instanton-induced four fermion interaction model in which the interaction is energy- 
independent; i.e., instantaneous. In addition, it should be noted that the above approach 
for getting the nonperturbative T 4 vertex has been proven to be very useful for the studies 
of nonperturbative vector and axial vector vertex [10- 11, 16]. This approach employs a con- 
sistent treatment of the dressed quark propagator G and the dressed vertex T 4 , which are 
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both determined from the effective quark-quark interaction by the rainbow DSE for G and 
the inhomogeneous ladder BSE for T 4 . 

With these two "phase" characterized by qualitatively different momentum-dependent 
quark propagator (24, 25), the GCM can be used to explore chiral symmetry restoration and 
phase transition between "Wigner" and "Nambu-Goldstone" phase. 

To explore the possibility of a phase transition one must consider the relative stability 
of the confined and deconfined phase by computing the /x-dependence of vacuum pressure 
difference(or "bag constant" [3]). It is equivalent to calculating the difference between the 
tree- level auxiliary-field effective action [17] evaluated with the "Wigner" solution charac- 
terised by B{p 2 ) = 0, and the "Nambu-Goldstone" solution characterised By B(p 2 ) ^ 0[12]: 



bm = p[g { NG) ] - p[gP\ 

A 2 f + C 2 + B 2 



4iV c 

3 
2^3 



r d A P f 
J sds J sin 2 xdx I In 



+ Re 



A' 2 f + C' 2 

sA 2 sin 2 X + C 2 + B 2 



f{°A ~ v'a) + (P4 + ¥)(<TC ~ o'c) 



sA' 2 sin 2 x + C' 2 



(26) 



+Re s sin 2 x {cta - ct'a) + (s 1/2 cosx + if^j (o"c* - cr'c) } > 



with 



A = [A(s) + 2ti2s 1/2 cosx^^}, C = s^cosxiAis) + 2i^s l/2 cosx^- L ] + ifiA(s), 



ds 



ds 



A' = [A'(s) + 2tfis 1/2 cos X ^^-], C = s 1/2 cos X [A'(s) + 2t f is 1/2 cos X ^^-} + tfiA'(s), 



B = B(s) + 2ifis 1/2 cosx 



ds 
dB(s) 
ds 



A' 



C 



s sin 2 X A 2 + C 2 + B 2 ' 



s sin 2 X A' 2 + C' 2 ' 



a 



s sin 2 xA 2 + C 2 + B 2 ' a ° s sin 2 X A' 2 + C' 2 ' 
B{ji) > indicates the stability of the confined(Nambu-Goldstone) phase and hence the 
phase boundary is specified by £>(/i) = 0. B(fx)/B(0) is plotted in Figs. 1-2. The scale is 
calculated to be B(0) = 0.131 ~ 0.152 (GeV) 4 (see Tabled), which can be compared with the 
value (0.145 GeV) 4 commonly used in bag-like models of hadron[18]. It is positive when the 
Nambu-Goldstone phase is dynamically favoured; i.e., has the highest pressure and become 
negative when the Wigner pressure become larger. The critical chemical potential is the zero 
of B(jj): i.e., n c = 270 ~ 375 MeV. This abrupt switch from the Nambu-Goldstone to the 
Wigner phase signals a first order transition. It should be noted that our numerical results is 



only valid for small values of chemical potential. If it could be qualitatively extrapolated to 
the reasonably large /x c , then our conclusion about the first order phase transition is correct. 




0.0 0.1 0.2 0.3 0.4 0.5 



H (GeV) 

Fig.l. Ratio B(fi)/B(0) for gluon propagator Eq.(14). 
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set 1 (A=10~ 1 GeV 4 , \ 

X=1.83GeV) \ 

set2(A=10" 4 GeV 4 , \ 

X=1.02GeV) \ \ 
set3(A=10" 7 GeV 4 , \ 1 
X=0.83GeV) \ 

i , i . i \ 
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Fig. 2. Ratio B(fi)/B(0) for gluon propagator Eq.(15). 



Table. I. The bag constant for two different models gluon propagator 



9 2 




2xi e £ 

A 2 




-A-K 2 d f 


A[GeV 2 } 


x[GeV] 


B*(0) [GeV] 


A[GeV 4 ] x[GeV] 


B$(0) [GeV] 


0.200 


1.65 


0.136 


10- 1 1.83 


0.135 


0.020 


1.55 


0.139 


10~ 4 1.02 


0.136 


0.002 


1.45 


0.131 


10~ 7 0.83 


0.152 



The chiral quark condensate is proportional to the matrix trace of the chiral-limit dressed 
quark propagator. Using Eq.(24), we obtain the following expression valid in the domain of 
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confinement and DCSB: 

(0| : qq : |0) M = -tr DC {Re . (27) 

The calculated ratio (0| : qq : |0) M /(0| : qq : |0) can be seen from the Figs. 3-4. In Figs. 3-4, 
we see that (0| : qq : |0) M decreases with increasing /z, up to a critical value of /i c , when 
it drops continuously to zero(It should be noted that although our numerical result is only 
valid for low values of chemical potential, we may still expect that it qualitatively gives the 
reasonable large /i tendency of the two-quark condensate). The change of (0| : qq : |0) M with 
increasing /z in our calculation is consistent with results in Ref. [19-20]. However, this result 
is qualitatively different from that in Refs. [13,14]. In Refs.[13,14], the chiral order parameter 
increases with increasing chemical potential up to /i c . At /i c , it drops abruptly to zero. We 
want to stress that the conclusion of Refs. [12-14] is based on the fact that those authors 
adopt the following general form for the solution of Eq.(8): 

g-\p) = %1 . p A(p) + i l4 (p 4 + ifi)C(p) + B(p), (28) 

where p = (p,p 4 + i/i). Although not explicitly indicated, the solutions(A(p), B{p) and C{p)) 
are functions only of \p\ 2 and (p 4 + i/i) 2 [14]. Comparing Eq.(24) with Eq.(28), it is easily 
seen that our solution(24) can not be included in the general form(28) given by Refs. [12-14]. 

1.00 
0.95 
o 0.90 

-i— » 
CO 

0.85 



0.80 

0.0 0.1 0.2 0.3 0.4 

n (GeV) 

Fig. 3. Ratio (6| : qq : |6) M /(6| : qq : |0) for gluon propagator Eq.(12). 




set 1 (A=0.2GeV 2 ,x=1-65GeV) 
set 2 (A=0.02GeV 2 , z =1.55GeV) 
set 3 (A=0.002GeV 2 , x =1 .45GeV) 

_i I i I i L 
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1.00 
0.95 
O 0.90 

-t— ■ 

2 

0.85 



0.80 

0.0 0.1 0.2 0.3 0.4 

n (GeV) 

Fig. 4. Ratio (0| : qq : |6) M /(6| : qq : |0) for gluon propagator Eq.(13). 

To summarize: in the present paper, we provide a general recipe to calculate the low 
chemical potential dependence of the dressed quark propagator at the mean field level in the 
framework of GCM. From this the /x dependence of the quark condensate and bag constant 
are evaluated. It is found that the chiral order parameter decreases with increasing /i, 
which is qualitatively different from that in Refs. [13,14]. The reason for this discrepancy is 
discussed. 
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